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polyhedra, the diamonds and the prisms. Except for the fullerene, all the spectra are 
obtained in closed form. They exhibit a rich pattern of degeneracies. The total energy 
at half filling is also evaluated in all the examples as a function of the magnetic charge. 
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1. Introduction 

Electronic properties of low- dimensional systems have been the focus of intense research 
for several decades [H [21 [3]. Electrons can be confined in various geometries, such 
as a plane (two-dimensional electron gas), a one- dimensional (quantum) wire, or a 
quantum dot, which can be referred to as a zero- dimensional system. In each one of these 
geometries, new and spectacular physical phenomena have been exposed: quantum Hall 
effect (in two dimensions), Luttinger liquid (in one dimension). Coulomb blockade and 
Kondo effect (in zero dimension), to mention just a few examples. 

The present work is concerned with another realization of a low-dimensional system 
which has so far received much less attention, namely where electrons are confined to 
move on a compact surface. The simplest class of such surfaces has the topology of 
the sphere. The fullerene and its derivatives provide the most natural candidates for 
such systems. The main goal of the present work is to investigate the single-particle 
energy spectrum as a function of a radial magnetic field, whereas the effect of spin-orbit 
interactions in the presence of a radial electric field is the subject of the companion 
paper [4]. 

Although the design of a pertinent experiment seems to be a formidable task, 
there are several motivations for pursuing this topic on the theoretical front. First, one 
encounters here a zero-dimensional system with a topology different from the usual one 
of a quantum box, say. It is surprising that, although the Cqo fullerene was discovered 
more than two decades ago [5] , the study of electronic properties of nanoscopic systems 
with the topology of the sphere is still rather scarce. It is indeed expected that topology 
and geometry should play an important role in the physics of electrons residing on 
closed surfaces. The examples treated in this work will demonstrate that the tight- 
binding spectra of polyhedra exhibit a rich pattern of degeneracies, at variance with 
the regularity of the spectrum of the continuum Schrodinger equation on the sphere 
in the presence of a magnetic charge [6]. A similar scenario was already known to 
occur in the case of the plane, when the simple pattern of Landau levels for an electron 
subject to a perpendicular magnetic field is compared with the strikingly rich Hofstadter 
butterfly [7] resulting when the same problem is treated on the square lattice, within 
the tight-binding scheme. Among more direct motivations, let us recall that it has 
been suggested that fuUerenes, i.e., carbon molecules with spherical topology, can be 
described in an effective way by the Dirac equation in the continuum in the presence 
of the magnetic field of an effective quantized monopole of charge one-half [8], in the 
absence of any physical magnetic field. Moreover, it was shown recently that magnetic 
monopoles may emerge in a class of exotic magnets known collectively as spin ice |9j. 
Third, the idea of considering electrons confined to a spherical shell in the continuum 
and subject to a magnetic field emerging from a magnetic monopole proved to be useful 
in the study of the quantum Hall effect [lOl [11], as it avoids the boundary problem. 

Having laid down our motivations, let us now move on and detail our work. The 
objective is to investigate the dependence of the tight-binding energy spectrum on the 
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value of a radial magnetic field, for a wide range of examples of graphs with spherical 
topology. More concretely, the electron lives on the vertices (sites) of a polyhedron 
drawn on the unit sphere and executes nearest-neighbor hopping. The energy levels 
are examined as a function of a radial isotropic magnetic field. The crucial difference 
between the present study and that of the Hofstadter problem [7] is that here, the 
only way to create a radial magnetic field is by placing a magnetic charge at the 
origin, i.e., at the center of the sphere. As is well known, such a magnetic charge 
must be quantized [12], i.e., be equal to the product of the elementary magnetic charge 
g = he/ (2e) of the famous Dirac monopole by an integer n. The total magnetic flux 
coming out through the sphere is n$0; where $o = 47r(7 = hc/e is the flux quantum. The 
(positive or negative) integer n, hereafter referred to as the magnetic charge, determines 
the gauge sector of the problem. Within the tight-binding model, the magnetic field 
enters the problem through U(l) phase factors living on the oriented links. The product 
of the phase factors living on the anticlockwise oriented links around a face of the 
polyhedron equals exp(27ri</9/$o)5 where ip is the outgoing magnetic flux through that 
face. This discrete formalism allows one to avoid the singularities met in the continuum 
description of monopoles [I3j. For a given graph drawn on the unit sphere, e.g. a regular 
polyhedron, the energy spectrum and all the observables depend solely on the magnetic 
charge n. Although it might be tempting to regard the value of n as representing the 
strength of the magnetic field, this viewpoint might be somewhat misleading. Recall that 
the Dirac quantization condition implies that there is no continuous deformation of one 
gauge sector into another. Furthermore, the results of this work will demonstrate that 
there is in general no simple relationship between the value of n and the observed pattern 
of level degeneracies. As already mentioned, the present situation is in strong contrast 
with its continuum analogue, namely that described by the Schrodinger equation on the 
sphere in the presence of a magnetic monopole, investigated in the pioneering work of 
Tamm [6], whose spectrum has a very regular dependence on the magnetic charge n. 
The energy levels indeed read 



with M being the mass of the particle, R the sphere radius, and £ = 0,1,... the angular 
quantum number. The corresponding multiplicities are 



The setup of the present paper is the following. The general concepts and notations 
are introduced in Section [2l In Section [3] the five regular polyhedra or Platonic solids 
(tetrahedron, cube, octahedron, dodecahedron, icosahedron) are considered. The F 
faces of these polyhedra are equivalent, so that the problem is periodic in the integer n, 
with period F. Section H] is devoted to the investigation of the Ceo fullerene, modeled as 
a symmetric truncated icosahedron (where all links have the same length) made of 12 
pentagons (with solid angle fl^) and 20 hexagons (with solid angle Og). The ratio fl^/fle 
is irrational, so that the spectrum depends on the magnetic charge n in a quasiperiodic 





nin/ = 2i + \n\ + 1. 
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fashion. Finally in Section two families of polyhedra are considered: the A^-prism 
(rectangular prism with a regular A^-gonal basis) and its dual the A^-diamond. The 
spectra of these two examples are explicitly derived for all > 3. Section [6] contains a 
discussion. Miscellaneous results on the geometrical characteristics of the Cgo fuUerene 
are presented in a self-consistent fashion in the Appendix. 



2. Generalities 



Throughout the following, the objects to be considered are polyhedral graphs drawn on 
the unit sphere. At the origin (center of the sphere) sits a quantized magnetic charge 
equal to n times the magnetic charge 5^ of a Dirac monopole, so that the total outgoing 
magnetic flux through the sphere is n times the flux quantum $o- 

Let us denote by V the number of vertices (sites), by L the number of links (bonds) 
and by F the number of faces of the polyhedron. In the case of the spherical topology, 
the Euler relation reads (see e.g. [H]) 

V-L + F = 2. (2.1) 

The tight-binding Hamiltonian of the system is 

Tt; = ^ (a\Uijaj + h.c.) , (2.2) 

<i3> 

where the sum runs over the L oriented links < ij > of the polyhedron, and the Uij 
are phase factors (elements of the gauge group U(l)) living on these links. The 
Hamiltonian is complex, i.e., not invariant under time reversal. One has 

= Uf^ = U% (2.3) 

where the star denotes complex conjugation. The product of the phase factors around 
each face is related to the outgoing magnetic flux ip through the face as 

UijUjk ■■■U^i = exp(27ri(^/$o). (2.4) 

The Hamiltonian Til is thus represented hj a V x V complex Hermitian matrix Ti 
such that Hij = Uij. The equation for the energy eigenvalues Ea and corresponding 
eigenf unctions 4'i,a reads 

Eai^i,a = (2.5) 

m 

with a = 1, . . . ,V, whereas j{i) runs over the neighbors of i. The energy spectrum 
obeys the following sum rules 

J2E, = 0, Y.El = 2L, (2.6) 

a a 

where the energy eigenvalues are counted with their multiplicities. The first sum 
indeed equals J^a^a = tiH = J^i'Hu = 0. This zero sum rule holds whenever 
the Hamiltonian only has non-diagonal matrix elements. The second sum equals 
J2aEl = tiTi? = I'^ijP = 2L, as each link gives two contributions equal to unity. 
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The main step in the exphcit construction of the Hamiltonian matrix Ti consists in 
finding a configuration of the gauge field, i.e., a set of phase factors f/jj, so that (12.41) 
holds for all the F faces. The L phase factors Uij are therefore constrained by F 
conditions. It will be clear from the example of the tetrahedron that only F — 1 of 
them are independent. Using (12.11) . one is then left with L — [F — 1) = V — 1 degrees 
of freedom in the choice of a gauge field. It proves convenient to fix these degrees of 
freedom by setting all the phase factors of the links of a spanning tree to their trivial 
values Uij = 1, along a line of thought dating back to [15]. It should be recalled that 
a spanning tree is a tree, i.e., a graph with no loop, that spans the graph, i.e., passes 
through every vertex of the graph (see e.g. [HI). The choice of a gauge thus boils down 
to the discrete choice of a spanning tree. The number of non-trivial phase factors is 
equal to its minimal value F — 1. 

As already announced, we shall successively investigate in detail the following 
examples: the five regular polyhedra or Platonic solids, the Cgo fuUerene (modeled as 
a symmetric truncated icosahedron) in Section HI and finally two families of polyhedra, 
the diamonds and the prisms, in Section [51 



3. Platonic solids 



There are only five regular polyhedra in three dimensions, the Platonic solids. This fact 
known since Antiquity can be easily recovered as follows. Let p be the coordination 
number of a vertex and q be the number of sides of a face. Evidently, the equality 
2L = pV = qF holds. The Euler relation (12. ip then implies 



1111 

- + - = o + 7- 
p q 2 L 



(3.1) 



The five solutions to (13.11) with p > 3 and q > 3 correspond to the Platonic solids, as 
recalled in Table [H The tetrahedron is its own dual, whereas (cube, octahedron) and 
(dodecahedron, icosahedron) form dual pairs (see e.g. |14j). 



polyhedron 


p 




V 


L 


F 


tetrahedron 


3 


3 


4 


6 


4 


cube 


3 


4 


8 


12 


6 


octahedron 


4 


3 


6 


12 


8 


dodecahedron 


3 


5 


20 


30 


12 


icosahedron 


5 


3 


12 


30 


20 



Table 1. The five regular polyhedra or Platonic solids, with their coordination 
numbers p, numbers q of sides of a face, and total numbers V of vertices, L of links 
and F of faces. 



The F faces of a regular polyhedron are equivalent. In particular, they have equal 
solid angles. The magnetic fiux through each phase is therefore = u^q/F, so that the 
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(3.2) 



(3.3) 



This property ensures that the problem is periodic in the integer n, with period F. 
3.1. The tetrahedron 

The tetrahedron, shown in Figure [1], is the simplest of the Platonic solids. It consists 
of 4 trivalent vertices, 6 links and 4 triangular faces. 



Throughout the following it will be advantageous to unwrap the polyhedra around 
an axis of high symmetry [16]. For the tetrahedron it is convenient to choose the 
threefold symmetry axis going through A and the center of the opposite BCD face. The 
planar representation thus obtained is shown in Figure [2l Each face appears exactly 
once, whereas vertices and links may have several occurrences, to be identified by the 
inverse procedure of wrapping the planar representation onto the sphere. 




A 



C 



Figure 1. The tetrahedron. The dashed-dotted Une shows the threefold axis used to 
unwrap the structure. 



A 



A 



• 




D 



Figure 2. Planar representation of the tetrahedron. Links drawn as dashed lines 
form the chosen spanning tree, whereas those drawn as full lines carry non-trivial 
phase factors. 



In order to construct the gauge field along the lines of the argument given in 
Section [2], we have chosen a convenient, highly symmetric spanning tree, shown in 
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Figure [2] as dashed lines. The phase factors of these hnks are set to their trivial values 
Uij = 1. The non-trivial phase factors are then obtained as follows. The condition (12.41) 
for the faces ABC, ACD and ADB respectively yields Ubc = ^, Ucb = ^ and Ubb = ^■ 
As far as the remaining face BCD is concerned, using (12.31) . one arrives at the relation 
UcbUbdUdc = The condition (12. 4p for the face BCD is automatically fulfilled, as 
a consequence of (13.31) which reads u;'^ = 1 since F = 4, so that uj~^ = uj. This explicit 
example illustrates the general properties that only F — 1 of the F conditions (12.41) are 
independent. To sum up, the trivial and non-trivial phase factors of the gauge field 
respectively read 

Uab = Uac = Uad = 1, f^Bc = UcD = Ubb = ^- (3.4) 
The Hamiltonian matrix thus obtained, 



n 



(0 


1 


1 


1 \ 


1 





UJ 




1 







UJ 




UJ 




J 



(3.5) 



depends on the magnetic charge n through the complex number uj introduced in (13. 2p . 
The energy eigenvalues of the above 4x4 matrix are listed in Table [2] and shown in 
Figure [3] for each value of the magnetic charge over one period. 



n 


E 


m 




3 


1 











-1 


3 


1,3 


v/3 
-73 


2 
2 




1 


3 


2 








-3 


1 



Table 2. Energy levels E of the tetrahedron and their multiplicities m, for each value 
of the magnetic charge n over one period. 



The other Platonic solids will now successively be dealt with along the same lines. 



3.2. The cube 

Figure H] shows the planar representation obtained by unwrapping the cube around a 
fourfold axis going through the centers of the opposite faces ABCD and EFGH, together 
with the spanning tree chosen to fix the gauge. One has uj^ = 1. The non-trivial phase 
factors of the gauge field read 

^FB = Ur,A = Ueu = UJ, Ugc = UuD = uj^- (3.6) 

The corresponding 8x8 Hamiltonian matrix is too complex to be diagonalizable 
by hand. For the cube and the subsequent Platonic solids, this task has been performed 
with the help of the software MACSYMA. The energy eigenvalues thus obtained and 
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12 3 4 



n 

Figure 3. Energy levels E of the tetrahedron as a function of the magnetic charge n 
over one period. Line segments show how individual levels 'jump' as n is increased by 
one unit. 




Figure 4. Planar representation of the cube. Same convention as in Figure [2l 



their multiplicities are listed in Table [3l and shown in Figure [5], for each value of the 
magnetic charge over one period. 



n 


E 


m 




3 


1 




1 


3 









-1 


3 




-3 


1 






2 


1,5 





4 






2 



n 


E 


m 




2 


3 


2,4 





2 




-2 


3 




V3 


4 


3 




4 



Table 3. Energy levels E of the cube and their multiplicities m, for each value of the 
magnetic charge n over one period. 
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1 2 3 4 5 6 
n 

Figure 5. Energy levels E of the cube as a function of the magnetic charge n over 
one period. Same convention as in Figure [3l 

3.3. The octahedron 

Figure [6] shows the planar representation obtained by unwrapping the octahedron around 
a fourfold axis going through the opposite vertices A and F, together with the spanning 
tree chosen to fix the gauge. One has uj^ = 1. The non-trivial phase factors of the gauge 
field read 

f^BC = UcD = Ube = Ueb = Uef = Ufc = t^^, f^FD = i^^- (3.7) 




F F F F 



Figure 6. Planar representation of the octahedron. Same convention as in Figure [2l 

The energy eigenvalues of the corresponding 6x6 Hamiltonian matrix and their 
multiplicities are listed in Table HI and shown in Figure [71 for each value of the magnetic 
charge over one period. 

3.4. The dodecahedron 

Figure [H] shows the planar representation obtained by unwrapping the dodecahedron 
around a fivefold axis going through the centers of the opposite faces ABODE and 
PQRST, together with the spanning tree chosen to fix the gauge. One has u^"^ = 1. The 
non-trivial phase factors of the gauge field read 

^EA = f^GB = f^QL = f^PT = ^, ^HC = ^RM = t^^, 

Uj^ = Us^ = iu^ Uje = Uto = uj\ U,K = io^ (3.8) 
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n 


E 


m 




4 


1 








3 




-2 


2 


1,7 


2V2 
-V2 


2 
4 




2 


3 


2,6 






-2 


3 



n 




m 


3,5 


V2 
-2v^ 


4 
2 




2 


2 


4 





3 




-4 


1 



Table 4. Energy levels E of the octahedron and their multiplicities m, for each value 
of the magnetic charge n over one period. 




Figure 7. Energy levels E of the octahedron as a function of the magnetic charge n 
over one period. Same convention as in Figure [S] 



D 



C J 

H .J ,.J 

I 
I 

S T *P 

I I 
I I 

R 



Figure 8. Planar representation of the dodecahedron. Same convention as in FigureO 



The energy eigenvalues of the corresponding 20 x 20 Hamiltonian matrix and their 
multiphcities are hsted in Table El and shown in Figure [HI for each value of the magnetic 
charge over one period. 
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n 


E 


m 




o 
o 


i 




/F 

V5 


3 




1 


5 








4 




-2 


4 






3 




V3(l + v5)/2 


z 




AT 

V3 


4 


1,11 


(v7- v3)/2 


n 





V3(l - V5)/2 


2 




-(n/7 + V3)/2 


6 




(3 + V5)/2 


3 




(V13- l)/2 


5 


2,10 


(3 - V5)/2 


3 




-1 


4 




-(v^+l)/2 


5 




V6 


4 




1 


6 


3,9 


-1 


6 






4 



n 


77" 


m 




(V13 + l)/2 


5 




1 


4 


4,8 


(V5- 3)/2 


3 




(1- v^)/2 


5 




-(3 + v^)/2 


3 




(V7 + V3)/2 


6 




V3(V5- l)/2 


2 


5,7 


(v^- V7)/2 


6 




/ — 


4 




-V3(l + v5)/2 


2 




Vo 


o 
o 




2 


4 







4 


6 


-1 


5 




-v^ 


3 




-3 


1 



Table 5. Energy levels E of the dodecahedron and their multiplicities m, for each 
value of the magnetic charge n over one period. 




10 12 



Figure 9. Energy levels E of the dodecahedron as a function of the magnetic charge 
n over one period. Same convention as in Figure [31 



3.5. The icosahedron 



Figure [TO] shows the planar representation obtained by unwrapping the icosahedron 
around a fivefold axis going through the opposite vertices A and L, together with the 
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spanning tree chosen to fix the gauge. One has uo"^^ = 1. The non-trivial phase factors 
of the gauge field read 



f^BC — f^CD — f^DE — UeF — UfB — 

Ugk = Uk3 = U,n = f/ffl = f^HG = ^, 

Ubg = Unc = ^ic = f^FG = Um = Ufk = 




Figure 10. Planar representation of the icosahedron. Same convention as in FigureO 



The energy eigenvalues of the corresponding 12 x 12 Hamiltonian matrix and their 
multiplicities are listed in Table El and shown in Figure [TT], for each value of the magnetic 
charge over one period. 




5 10 15 20 



n 

Figure 11. Energy levels E of the icosahedron as a function of the magnetic charge 
n over one period. Same convention as in Figure [S] 



3. 6. Properties of the spectra 

The energy spectra of the five Platonic solids have been given in Tables |2] to [6] and 
illustrated in Figures [3l [5l [TJ [9] and [11], as a function of the magnetic charge n over one 
period. These energy spectra exhibit several remarkable properties. First of all, the 
energy levels are even functions of the magnetic charge n, i.e., they are invariant under 
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(J 


75 
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-1 


5 




-v^ 


3 




a/10(5 + a/5)/2 


2 


1,19 


\lb - 2^/5 


4 




-\ 2ih + \/5)/2 


6 


2,18 


(5 + V5)/2 
(v^-3)/2 


3 
5 
4 




\/5 + 2\/5 


4 


3,17 




6 




-^10(5- v^)/2 


2 


4,16 


(3 + V5)/2 
(v^-5)/2 
-v^ 


5 
3 
4 


5,15 


v/5 


6 
6 



n 


E 


m 




V5 


4 


6,14 


(5 - Vh)/2 


3 




-(3 + v^)/2 


5 




x/l0(5- v^)/2 


2 


7,13 


a/2 (5 - v^)/2 


6 




+ 2\/5 


4 




v/5 


4 


8,12 


(3 - V5)/2 


5 




— (,0 + V oj/^ 


Q 




^^2(5 + V5)/2 


6 


9,11 


— Y "3 — V 3 


4 




-^/l0(5 + v^)/2 


2 




V5 


3 


10 


1 


5 


-v^ 


3 




-5 


1 



Table 6. Energy levels E of the icosahedron and their multiplicities m, for each value 
of the magnetic charge n over one period. 



a change of the sign of the magnetic field. As announced above, they are also periodic 
in n, with period F . 

Another symmetry of the energy spectra can be revealed as follows. Suppose all 
the (trivial and non-trivial) phase factors Uij are simultaneously changed into their 
opposites. The whole Hamiltonian matrix Ti, and therefore all the energy levels, are 
also changed into their opposites. On the other hand, as all the faces have q sides, the 
product in the left-hand side of (12. 4p is multiplied by (—1)'. The magnetic flux through 
each face has therefore been increased by 5(p = g$o/2, i.e., n has been increased by 
qF/2. Therefore, if q is even (this is the case for the cube), the energy spectrum is its 
own opposite, i.e., it is symmetric with respect to the origin of energies, for every value 
of n. If q is odd (this is the case for the four other Platonic solids), the spectrum at 
n + F/2 is the opposite of that at n, i.e., the spectrum is a semi-periodic function of n. 

All the energy levels have been found to be algebraic numbers of degree at most four, 
with relatively simple expressions involving square roots. This property is somewhat 
surprising, especially in view of the dimension of the Hamiltonian matrices, which can 
be as high as \^ = 20 for the dodecahedron. The energy levels fulfill the sum rules (12. 6p . 
Checking this expected property is however a non-trivial exercise in algebra for the case 
of the icosahedron. 
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The multiplicities of the energy levels are typically rather high, and they vary 
in a complex and irregular fashion as a function of both the magnetic charge and 
the rank of the considered energy level. The spectra however exhibit the following 
property: the largest energy eigenvalue has multiplicity m = |n| + 1 for a small enough 
absolute magnetic charge < rimax, where rimax = 2,3,3,5,5 for the five Platonic 
solids ordered as above). This observed linear increase is a perfect analogue of the 
continuum result (11.21) . with the reasonable assumption that the continuum counterpart 
of the largest energy eigenvalue corresponds to the orbital quantum number £ = 0. 



3.7. Total energy 

An interesting illustration of the above energy spectra is provided by the total energy 
at half filling, defined as 

V/2 

S = T.Ea, (3.10) 

a=l 

where the V energy levels are assumed to be in increasing order {Ei < E2 < ■ . . < Ey) 
and to be again repeated according to their multiplicities. The first of the sum 
rules (12.61) implies that the total energy £ thus defined is insensitive to the sign of 
the Hamiltonian 7i, i.e., it would be unchanged if the matrix elements were replaced by 
their opposites Tiij = —Uij. The second of the sum rules (12.61) implies that the mean 
squared value of the individual energy levels reads {E^) = 2L/V = p, where p is the 
coordination number of the vertices. This observation is indicative that the total energy 
scales as y/pV = and therefore suggests to introduce the reduced total energy 

£ £ 

£r=-^ = ^=. (3.11) 

This heuristic argument can be turned to a quantitative prediction in the limit of 
an infinitely connected structure (p — cx)), such as e.g. the hypercubic lattice in high 
spatial dimension d, where p = 2d. Although the full energy spectrum is the interval 
—p ^ E < p, typical energy eigenvalues only grow as ~ y/p in the limit under 
consideration, and the bulk of the normalized density of states becomes Gaussian: 

The reduced total energy therefore has the following universal limiting value: 

1 /"O 



i.e., 



Ep^{E)dE, (3.13) 



£00 = i= = -0.398942. (3.14) 

V 27r 

Figure [12] shows a plot of the reduced total energy £,. as a function of the reduced 
fiux per face, n/F = f/^o, for the five Platonic solids. This quantity is observed to 
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vary within a rather modest range around the hmiting value fl3.14l) . shown as a dashed 
hne. 

In the case of the cube, is minimal for n/F = 1/2, in agreement with the 
prediction that the total energy is minimal when the flux per face equals the filling 
factor [17] . In the four other cases, however, the last of the above symmetry properties 
implies that Sj- is periodic in n with period F/2, and not F. As a consequence, it takes 
the same value for n/F = 1/2 and in the absence of a magnetic field. As it turns out, Sj- 
is observed to be minimal for n/F = 1/4 and 3/4 (F is a multiple of 4). 

-0.25 1 ' ' ' ' ' ' ' 1 



-0.35: 



-0.45 




■ 0.25 0.50 0.75 1.00 

n/F 

Figure 12. Plot of the reduced total energy £r against the reduced flux per face n/F, 
for the five Platonic sohds: tetrahedron (down triangles), cube (stars), octahedron (up 
triangles), dodecahedron (squares) and icosahedron (circles). The dashed line shows 
the hmiting value p.l4p . 



4. The Ceo fullerene 

We now turn to the case of the Ceo fullerene. This carbon molecule with the shape 
of a truncated icosahedron has been discovered in 1985 For simplicity we model 
it as a symmetric truncated icosahedron, where all the links have equal lengths. This 
symmetry is known to be slightly violated [18], as for the free molecule the length of 
the sides of the pentagons is 1.46 A, whereas the length of the other links is 1.40 A. 
Considering the symmetric polyhedron will however not affect the salient qualitative 
features of our results. It will indeed turn out that these features can be explained by 
the quasiperiodic dependence of the energy spectrum on the magnetic charge. 

The symmetric truncated icosahedron has = 60 equivalent vertices, L = 90 
equivalent links, and F = 32 faces, namely 12 pentagons and 20 hexagons, respectively 
corresponding to the vertices and to the faces of the icosahedron. Various results on 
its geometrical characteristics, which are useful both in the present work and in [1], 
are presented in a self-consistent fashion in the Appendix. Figure [13] shows the planar 
representation obtained by unwrapping the fullerene around a fivefold axis going through 



Tight-binding spectra on spherical graphs I: the effect of a magnetic charge 



16 



the opposite pentagonal faces A1A2A3A4A5 and L1L2L3L4L5, together with the spanning 
tree chosen to fix the gauge. The labehng of the vertices is consistent with that used in 
Figure [THl 




Figure 13. Planar representation of the fuUerene. Same convention as in Figure [2l 
The 12 pentagonal faces Ai . . .A5 to Li . . .L5 have not been drawn, for clarity. 



The magnetic fiuxes ^p^ and y^e through a pentagonal and a hexagonal face are 
proportional to the solid angles VL^ and fig, given by (lA.14p . If n denotes the magnetic 
charge, the corresponding phase factors read 

a = exp(27ri(y95/$o) = exp(inr25/2), 

P = exp(27riv96/$o) = exp{inQ(i/2). (4.1) 
Equation flA.12|) implies = 1. The non-trivial phase factors of the gauge field are 





= f^LiLs 


= UbsB4 


= f^G4G3 


= f^C4C3 = ^D4D3 = OC, 






UE4E3 


= ^F4F3 


= ^H4H3 


= f^l4l3 = 


= f^J4J3 = f^K4K3 = Oij 








= UC2D5 


= UD2E5 


= UE2F5 




- ^HaGs 


= P. 


UF2B5 


= f^GaKs 


= af3, 


%4G4 = 


t/H4C4 = 








= a'(3\ 


Ul4Y)4 = 




t^J3D3 = a'^P^, Uj^E4 = 






^KgEs 


= a'p'' 


t^K4F4 


= a'p'' 


, [/g3F3=«V. 




(4. 



Figure [Ml shows a plot of the energy spectrum of the fullerene against the magnetic 
charge n, up to n = 300, obtained by means of a numerical diagonalization of the 
Hamiltonian matrix. In the absence of a magnetic monopole, we recover the known tight- 
binding spectrum of the fullerene [19], with its 15 different levels having multiplicities 
ranging from 1 to 9. For a non-zero magnetic charge n, the observed pattern of 
multiplicities only depends on the parity of n: 

• If n is even, the spectrum consists of 16 distinct energy levels: 1 nondegenerate (with 
multiplicity 1), 6 with multiplicity 3, 4 with multiplicity 4 and 5 with multiplicity 5. 

• If rz is odd, the spectrum consists of 14 distinct energy levels: 4 with multiplicity 2, 
4 with multiplicity 4 and 6 with multiplicity 6. All the multiplicities are even in 
this case. 
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Finally, the multiplicity of the largest energy eigenvalue grows linearly with the 
absolute magnetic charge, according to the rule m = |?t,| + 1, already observed in 
the case of the Platonic solids, up to nmax = 5 in the present case. 





Figure 14. Energy spectrum of the fuUerene as a function of the magnetic charge n. 
The arrows show the peaks at multiples of the quasi-period Q = 84. 



The Hamiltonian 7^, whose matrix elements involve a and /3, is a 47r-periodic 
function of nVt^ and nfig, viewed as two independent variables. In other words, it 
is a quasiperiodic function of the magnetic charge n, because the solid angle ratio 
R = Q^/Qq, given by (1A.16I) . is an irrational number. The quasi-period Q = 84, 
shown as arrows in Figure [TH clearly emerges from the pattern, especially as sharp 
peaks where the highest energy level is very close to its value E = p = 3 in the absence 
of a magnetic charge. The quasi-period Q can be identified as the smallest integer such 
that the reduced magnetic fluxes through each type of face are very close to integers: 

^ = 9^ = 1.972412 ^ 2, 
$0 47r 

^ = ^ = 3.016552 ^ 3. (4.3) 

$0 47r ^ ^ 

The ratio 2/3 is the first of the sequence of best rational approximants to the solid angle 
ratio R given in (]A.16p . The next approximant, 17/26, corresponds to the quasi-period 
Q = 724. 

Figure [15] shows a histogram plot of the density of states obtained by accumulating 
the spectra up to = 10^. This procedure amounts to performing a uniform averaging 
over the independent angles nfl^ and uQq. First of all, the density of states is observed 
to be an even function of the energy E. This property can be easily explained along the 
lines of Section [321 The Hamiltonian H. is indeed a semi-periodic function of nfl^. 
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The very irregular behavior of the density of states, with its many narrow peaks 
whose height keeps growing as the bin size 5E is decreased, is reminiscent of the singular 
spectra of more conventional quasiperiodic structures, especially in low dimension, such 
as the Fibonacci chain (see [20j for a recent review). It is, however, worth noticing 
that the density of states also has a smooth and rather uniform background all over the 
spectrum, i.e., for —?>< E <?>. 

An intriguing question is whether the study of the fullerene spectrum can tell us 
something about the spectrum of graphene [H]. In the absence of magnetic field, the 
density of states of graphene vanishes linearly for energies close to half filling [21j, as 
p{E) oc \E\, whereas in the presence of a magnetic field B there are Landau levels 
at Ejn = sign(m)y2e/w^5]m[ (m = 0, ±1, ±2, . . .), where v is the band velocity j22j . 
In particular, there is a Landau level at zero energy, as Eq = 0. It is, however, too 
speculative to interpret the peak ai E = visible in Figure [T5] as a zero-energy Landau 
level. 




-3-2-10 1 2 3 
E 

Figure 15. Histogram plot of the density of states p{E) of the fullerene. Spectra up 
to a magnetic charge n — 10^ are accumulated. The bin size is SE = 0.02. 

Figure [TB] shows a plot of the reduced total energy at half filling, = S/{60y/3) 
against the magnetic charge n, up to n = 300. The data exhibit oscillations at the 
above quasi-period Q = 84, with a strong third harmonic, especially near the minima. 
The mean reduced total energy, (S^) = —0.4394, is only 10.2 percent larger (in absolute 
value) than the limiting value E^o of (13.141) . although the density of states shown in 
Figure [15] is very far from being a Gaussian. 

5. The diamond and the prism 

We end up this investigation by considering the following two families of polyhedra. 

• The diamond is obtained by gluing together two pyramids whose basis is a regular 
polygon with sides. It has + 2 vertices, the two poles (with coordination 
number A^) and A^ vertices along the polygonal equator (with coordination number 
4), 3A^ links and 2A^ equivalent triangular faces. 
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50 100 150 200 250 300 

n 

Figure 16. Plot of the reduced total energy £,■ of the fuUerene against the magnetic 
charge n. The arrows show that the minima are separated by multiples of the quasi- 
period Q = 84. The dashed line shows the mean (£r) = —0.4394. 



• The prism is a rectangular prism whose basis is a regular polygon with sides. It 
has 2N equivalent vertices with coordination number 3, 3A^ links, and + 2 faces, 
the two polygonal bases and lateral rectangular faces. 

There is an infinite family of each kind of polyhedra, labeled by the integer > 3. 
The A^-prism and the A^-diamond are dual to each other (see e.g. [1_4]). Two of 
the Platonic solids, namely the octahedron and the cube, are recovered for A^ = 4. 
Figures [T71 and [T8l respectively show the planar representations obtained by unwrapping 
the diamond and the prism around their A^-fold axis, together with the spanning trees 
chosen to fix the gauge. The tight-binding spectra will be worked out explicitly for each 
family of graphs, for arbitrary A^ and an arbitrary magnetic charge. 





C C C C 

Figure 17. Planar representation of the diamond. Same convention as in Figure [2l 



5.1. The diamond 

The 2N triangular faces of the diamond are equivalent. As a consequence, if the 
magnetic charge is n, the magnetic flux per face is n$o/(2A^), so that the right-hand 
side of (12. 4p reads 

u = exp{i7Tn/N) (5-1) 
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A, 



^N-l 



A, 




Figure 18. Planar representation of the prism. Same convention as in Figure [D 
The upper and lower polygonal bases Ai . . .Ajv and Bi . . .B^r have not been drawn, 
for clarity. 



and obeys 

o;^ = (-1)'^, = 1. (5.2) 

The latter property ensures that the problem is periodic in the integer n, with period 2N. 
With the notations of Figure [13, the non-trivial phase factors of the gauge field 

read 

f^B,B,+, t/cB, =^'' (5.3) 

for k = 1,...,N, with periodic boundary conditions (A^ + 1 = 1). The eigenvalue 
equation (I2.5P therefore reads 

N N 

E^Pa = E ^k, E^c = E ^~^'^k, (5.4) 

fc=l k=l 

where the ipk = i^Bk obey periodic boundary conditions. Setting = u^'ipk, the above 
equations simplify to 

Eak = Ofc-i + Ofc+i + uj~''ipA + i^'^'ipc, 

N N 

E'ipA = E ^''^k, Eipc = E ^'''(^k, (5.5) 

k=l k=l 

where the obey the boundary conditions oq = iv'^aN, aN+i = cu^ai. The eigenvalues 
and eigenf unctions solving (15.51) can be derived exphcitly as follows. 

Consider first a generic value of the magnetic charge (n 7^ and n 7^ A^), so that u 
is not real, and look for an extended plane- wave solution of the form ak = e''^^, so that 
E = 2 cos q. The boundary conditions yield q = (2m — n)'n'/N, with m = 1, . . . , N. One 
has t/'a = V'c = 0, which imposes m ^ n and m ^ N . We thus obtain a band of — 2 
eigenvalues, 

(2m — n)7r 

E = 2cos^ (5.6) 



with a hole at 



Eo = 2cos^, (5.7) 
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corresponding to m = n being forbidden. Four other solutions correspond to impurity 



states: either ak = uj^, so that ipp^ = but ipc 0) '^fc 



UJ 



so that tpc = ^ but 



ipA 7^ 0. Both cases yield two minibands of twice degenerate impurity levels: 
E 



nn / nil 
cos— ±^iV + cos2— . 



(5. 



Consider now the values n = and n = N oi the magnetic charge, respectively 
corresponding to an integer and a half- integer flux quantum per face, i.e., to u = 1 and 
UJ = —1. There are now — 1 eigenvalues in the band (15.61) . which has no hole. Three 
other solutions correspond to impurity states. First, = and ^/'a + = yields one 
impurity state right at the band center, i.e., at energy 

E = 0. (5.9) 

Second, = and ipA = i'c yields two nondegenerate impurity levels: 

(5.10) 



E 



to ± y/2N + 1, 

lying further away from the band than the minibands of (15.81) . The energy spectrum 
thus obtained can be checked to obey the sum rules (12.61) with L = 3N, and to give 
back the spectrum of the octahedron (see Section 13. 3p in the case = 4. 

Figure [12] shows a plot of the energy levels of the diamond for = 30 as a function 
of the magnetic charge n over one period. The main features of the spectrum are clearly 
visible: the band (15. 6p of extended states with its hole ( 15. 7p . the minibands (15. Sp of 
impurity states, and the isolated impurity levels in the two commensurate cases, shown 
as larger dots, namely (15.91) within the band and (I5.10p away from the band. 

10 
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Figure 19. Energy levels E of the diamond for = 30 as a function of the magnetic 
charge n over one period. Large dots show the isolated impurity levels (|5.9p and (|5.10p . 
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5.2. The prism 

The prism has two different types of faces, two polygonal bases and rectangular 
faces. The solid angles of each type of face viewed from the center of the solid could 
be evaluated explicitly, using results from the Appendix, as a function of N and of the 
aspect ratio h/a, where h is the height of the prism and a the side of its polygonal bases. 
This parametrization will not be needed in the following. 

Let n be the magnetic charge, ^p-Q the flux through either of the polygonal bases, 
and if the flux through a rectangular face. One has then 2ip-Q + Nip = n$o- The 
corresponding phase factors, 

a = exp(27ri¥?B/$o), /? = exp(27ri¥?/$o), (5.11) 

obey 

a^/?^ = 1. (5.12) 

In the following we take /3 as our basic variable, at fixed magnetic charge n, forgetting 
about its dependence on the geometry of the prism. We introduce for further convenience 
the variable 

u = ^, (5.13) 

such that P = exp(i'u). 

With the notations of Figure [18], the non-trivial phase factors of the gauge field 

read 

t^AjvAi = %iBjv = «> f^B.A, = (5.14) 

The eigenvalue equation (12.51) therefore reads 

Eak = ttk-i + ttk+i + P''bk, Ebk = bk-i + bk+i+ P^'^ak, (5.15) 

where the = ipA^ ^i-nd bk = ipBk obey the boundary conditions 

bN = abo, b^+i = Cibi. Setting Ck = P'^bk, the above equations simplify to 

Eak = ak-i + ak+i + Ck, Eck = Pck-i + P~^Ck+i + ak, (5.16) 

whereas the Ck obey the same boundary conditions as the a^, i.e., Cq = aajy, Ci = aajy+i. 
The eigenvalues and eigenfunctions solving (15.161) can be derived explicitly as follows. 
Looking for an extended plane-wave solution of the form 

Equation f l5.16p implies that the energy E fulfills the condition 

2 cos q — E 1 

1 2 cos(g — u) — E 

Introducing the shifted momentum 

« = .-f = «-f. (5.19) 



0. (5.18) 
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the above equation yields the following two-band dispersion relation 

VL I IT 

E±{Q) = 2cosgcos- ± Wl + 4sin2Qsin2 -. (5.20) 

The boundary conditions lead to the following quantized values of the momentum: 

^ (2m + n)7r , , , , 

Q = - (m = l,...,iV). (5.21) 

Each energy level is an even and periodic function, whose period is, as expected, 2tt 
in u, i.e., $o in V^- Furthermore, the spectrum bears an extra discrete dependence on the 
magnetic charge n through the quantization condition (15.211) . In fact, only the parity 
of n matters, and levels for even n and odd n alternate. Figure [20] shows a plot of the 
energy spectrum of the prism with = 30 against u/{2'k) over one period, for an even 
magnetic charge n. The levels for odd n would lie between these levels. The energy 
spectrum thus obtained can be checked to obey the sum rules (12. 6p with L = 3N, and 
to give back the spectrum of the cube (see Section [3^ for = 4 and u = mi/3. 




0,2 0,4 0,6 0,8 1,0 

u/(2t\) 

Figure 20. Energy levels E of the prism for = 30 against u/ (27r) — (p/^o over one 
period, for an even magnetic charge n. 

The reduced total energy at half filling, defined in (13. lip , has a well-defined limit 
St.{u) as the number of sites of the prism becomes infinitely large, irrespective of the 
parity of the magnetic charge n. This limiting function has two different expressions, 
involving incomplete elliptic integrals, according to whether the bands E±{Q) overlap 
{u < 2n/3 and u > 47r/3) or not (27r/3 < u < 47r/3). Figure [21] shows a plot of S^{u) 
against u/{27i). The mean reduced total energy, (Sr) = —0.43866, is only 9.9 percent 
larger (in absolute value) than the limiting value S^o of (13.140 . 
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Figure 21. Plot of the limit £-[{u) of the reduced total energy of the prism, against 
u/{2it) = f/^o over one period. Dots show the values of u which separate the regions 
where the bands overlap or not, so that fr(u) has different analytic expressions. The 
dashed line shows the mean {£,■) — —0.43866. 



6. Discussion 

We have presented an extensive study of the tight-binding spectra on various polyhedral 
graphs drawn on the unit sphere, as a function of the magnetic field produced by a 
quantized magnetic charge sitting at the center of the sphere. For a fixed polyhedron, 
there is only one discrete parameter left, namely the integer magnetic charge n, fixing 
the gauge sector of the model. 

The spectra of the five Platonic solids, described in Section [3l exhibit a periodic 
dependence on the magnetic charge n, the period being the number of faces F of the 
polyhedron. The observed multiplicities of the energy levels are typically high. The 
multiplicity m = |n| + 1 of the largest energy eigenvalue, observed for the smallest 
values of the absolute magnetic charge, is a discrete analogue of the known continuum 
spectrum of multiplicities dating back to Tamm [6j. All the energy levels have rather 
simple expressions, involving at most two square roots, although the dimension of the 
Hamiltonian matrices can be as large as 20 in the case of the dodecahedron. It would 
be worthwhile to explore this rich pattern of multiplicities within a more formal group- 
theoretical framework. The case of the Cgo fullerene, modeled as a symmetric truncated 
icosahedron, is dealt with in Section HI The main features of the spectrum are explained 
by the quasiperiodic dependence of the energy spectrum on the magnetic charge n, 
and mainly by the occurrence of a quasi-period Q = 84, related to the first rational 
approximant of the area ratio. These features would still be present in the generic case 
of a non-symmetric truncated icosahedron, with two different bond lengths. Finally, the 
two families of polyhedra investigated in Section [5] have given to us the opportunity of 
underlining yet other features of the problem pertaining to less symmetric situations, 
such as the relevance of one-dimensional band structures in the very anisotropic regimes 
of the prism and the diamond at large N. 



Tight-binding spectra on spherical graphs I: the effect of a magnetic charge 25 

The investigation of physical properties has been focused onto the total energy £ 
at half filling. For all the polyhedra where all the vertices have the same coordination 
number p (i.e., all the examples considered in this work except for the diamond), the 
typical value of the total energy is found to be rather close (within a range of ten to 
twenty percent) to its universal large-p behavior corresponding to the asymptotically 
Gaussian nature of the density of states. The behavior of the total energy as a function 
of the magnetic charge n depends on the underlying polyhedron. The total energy is 
minimal when the magnetic flux per face equals 1/2 flux quantum for the cube, and more 
generally for the family of prisms, in agreement with the prediction that the minimum 
occurs when the flux per face equals the filling factor [17J. In the other four polyhedra, 
to the contrary, two energy minima are attained for the fractions 1/4 and 3/4, as a 
consequence of an extra symmetry of the spectra. 
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Appendix A. The Cgo fullerene. Geometrical characteristics 

For simplicity we model the Ceo fullerene as a symmetric truncated icosahedron, where 
all the links have equal lengths. This polyhedron consists of 12 pentagons and 20 
hexagons, respectively corresponding to the vertices and to the faces of the icosahedron. 
The goal of this Appendix is to derive in a self-contained fashion expressions for some 
of its geometrical characteristics, to be used both in this work and in |1] , and especially 
the solid angles Vt^ and VLq. 

A 




Figure Al. The triangular face ABC of the icosahedron decorated by vertices of the 
fullerene. Notations are consistent with Figures [TOl and fT3l 



Figure lAll shows an enlargement of the upper left part of Figures [TO] and [131 with 
consistent notations. We denote by A the vector from the origin to the vertex A, and so 
on. Choosing for convenience a coordinate system such that A is at the North pole and B 
lies in the x^-plane, the coordinates of the vertices A, B and C of the icosahedron read 

A =(0,0,1), 
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B = {sm9i, 0, cos6'i), 

C = (sm6'iCos(27r/5),sin6'iSm(27r/5),cos6'i), (A.l) 

Let 6i be the arc length of the hnks of the icosahedron. One has cos 6^ = A - B = B ■ C, 
hence 

cos^i = ^. (A.2) 
V5 

The vertices Ai, A2, Bi and Ci of the fuUerene can be parametrized as 
Ai = XB + jjA, A2 = AC + /iA, 

Bi = \A + fxB, Ci = \A + fxC. (A.3) 
The conditions Af = 1 and Ai ■ A2 = Ai ■ Bi yield n = 2X and 



_ 25-4V/5 

V 109 ^ ^ 

Let 6f be the arc length of the links of the fullerene. One has 

/7^^2 80 + 9^5 

COS = Ai ■ A2 = (4 + V5)A2 = — — — . (A.5) 

iuy 

The above expressions can be used to evaluate the relevant solid angles. Consider 
a spherical triangle defined by three points A, B, C on the unit sphere. It is well- 
known that the area (solid angle) of the triangle is given by the spherical excess 
Q = A + B + C — vr, where A is the angle of the triangle at its vertex A, and so 
on. This result is however not very useful in the present situation, where the vertices 
are known through their coordinates. A more convenient expression reads 



nfi s s — a s — b s — c /A^N 

tan — = tan -tan tan tan , A.o 

4 2 2 2 2 ^ ^ 

where a, b, c are the arc lengths of the sides, such that cos a = B ■ C , and so on, and 

s = (a + 6 + c)/2 is half the perimeter. The formula (1A.6I) is the spherical analogue of 

a well-known expression for the area of a planar triangle, 

= s{s-a){s-b){s-c), (A.7) 

dating back to Antiquity and known as Heron's formula. Equation (1A.6P can therefore 
be referred to as the spherical Heron formula. It can be recast, using trigonometric 
identities, into the following alternative form [24j : 

n I -I- cos^ ^ + cos^ I - 1 cosa + cos6 + cose + 1 , , 

cos — = ^ ^ ^ = r . (A. 8) 

2 aba aba 

2 cos - cos - cos - 4 cos - cos - cos - 

2 2 2 2 2 2 

Some further algebra led us to the following form, that is especially convenient in the 

case where the vertices are known through their Cartesian coordinates: 

"- (i + a.b)(i^;b^C)(h-c.A) - 
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where 

{A,B,C) = {Ax B) - C = A - {B X C) (A.IO) 

is the scalar triple product of the three unit vectors. 

For the isosceles triangle AA1A2, the formula (]A.9P yields after some algebra 



13 - 4v^ + V605 + 184v^ 

cosf2T = ^ . (A.ll) 

3d 

The solid angles of a pentagonal and a hexagonal face obey 

12^5 + 20r]6 = 47r (A.12) 

and read 

^5 = BQt, ^6 = - - 31^T- (A.13) 
5 

Some further algebra using trigonometric identities yields 

54887 + 720v^ ^ 511 + 945^5 

cosiZs = , cosiZfi = , A. 14 

^ 59049 ' ^ 2916 ' ^ ^ 

hence the numerical values 

^5 = 0.295072, ^6 = 0.451275. (A.15) 
To close up, it is worth comparing the ratio 

= ^ = 0.653863 (A.16) 

to its planar analogue, namely the area ratio of a pentagon and a hexagon with the 
same side. In the plane, the area of a regular n-gon of side a is 

An = , y (A.17) 

4tan(7r/nj 
The area ratio therefore reads 



A5 V15 5 + 2v^ ^ 
^ ^ 5 ^ V — ^ ^ 0.662212. A.18 

A 18 ^ ' 

The comparison of both results flA.16p and flA.18p shows that the effect of curvature is 
rather weak, as one has R/r = 0.987392. 
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